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A Langevin equation of heavy quarks in high-temperature quark-gluon plasma is derived. The 
dynamics of heavy quark color is coupled with the phase space dynamics and causes a macroscopic 
superposition state of heavy quark momentum. Decoherence of the superposition state allows one 
to use a classical description. The time scale of decoherence gives an appropriate discretization time 
scale At ~ for the classical Langevin equation, where M is heavy quark mass and 7 is 

heavy quark momentum diffusion constant. 
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I. INTRODUCTION 

The transport properties of quark-gluon plasma 
(QGP) have attracted a lot of interest since nearly per¬ 
fect liquid behavior has been discovered in the relativistic 
heavy-ion collision experiments at the Relativistic Heavy 
Ion Collider (RHIC) [1]. Further investigations in the 
heavy-ion collisions at higher energy are ongoing at the 
Large Hadron Collider (LHC). The discovery of nearly 
perfect liquid behavior seems to contradict the notion of 
weakly interacting QCP, and it rather suggests a strongly 
interacting nature of the QCP. Indeed, the universal 
lower bound of the ratio of shear viscosity to entropy 
density ??/s > l/47r is proposed in the strongly interacting 
quantum field theories based on the conjectured duality 
between the gauge theory and string theory, or the anti- 
de-Sitter space and conformal field theory (AdS/CFT) 
correspondence [2]. This is close to the value 77/5 ~ (1- 
3)/47r, which is phenomenologically extracted from the 
experimental data by hydrodynamic simulations for the 
heavy-ion collisions [3]. 

In the heavy-ion collisions, not only the bulk collective 
dynamics but also hard probes, such as heavy quarks and 
jets, reveal independent aspects of the transport proper¬ 
ties of the QCP. For example, medium modification of 
heavy quark momentum spectra offers an opportunity to 
study drag force of the QCP acting on a test particle 
with color [4]. Phenomenological studies of the heavy 
quark probe in the heavy-ion collisions are summarized 
in Ref. [5]. Since there is a kinematical hierarchy between 
the heavy quark and the QCP, the heavy quark dynamics 
is slow compared to the correlation time of matter con¬ 
stituents so that the effects of interaction between them 
are averaged out. This enables one to use a simple effec¬ 
tive description of the heavy quark using the drag and 
fluctuation forces. There have been various efforts to 
calculate the heavy quark transport coefficients by the 
perturbation theory [4, 6], by the lattice QCD simula¬ 
tions [7], and by applying the AdS/CFT correspondence 
[ 8 ]. 

In this paper, I shed light on the dynamics in the heavy 
quark color space. The Langevin dynamics is a classical 
description of the Brownian motion. Unlike kinetic vari¬ 


ables, the time scale of the heavy quark color does not 
get long even with its heavy mass. Typically, the time 
scale of color diffusion is about 1/g^T while the dura¬ 
tion of a soft scattering is about IjgT. Therefore, it is 
only in the weak coupling regime that one can expect 
the validity of Langevin description that couples with 
the heavy quark color degrees of freedom. In the case 
of the heavy quark with color, description in the phase 
space is classical while that in the color space is quan¬ 
tum. Here the heavy quark color is in the fundamental 
representation of color SU(Ac). Instead of describing it 
classically by Wong’s equation [9] in terms of a (-1)- 
dimensional vector, I treat a color state as a quantum 
state in a Hilbert space with dimensions. Indeed, the 
latter description is a direct consequence of the quantum 
description of the heavy quark Brownian motion [10]. 

The main finding of this paper is that an analog of 
the Schrodinger’s cat state appears in the heavy quark 
Langevin dynamics due to the non-Abelian interaction 
of QCD. The classical momentum corresponds to the cat 
and couples to quantum states in the color space. It is 
well known that decoherence of such a macroscopic su¬ 
perposition state is essential to have a classical picture of 
the Brownian motion [11]. Consequently, the discretiza¬ 
tion time scale of the classical Langevin dynamics must 
come out from the time scale of decoherence in the heavy 
quark sector. The discretization time scale turns out to 
be At ~ It depends on the heavy quark mass 

M and on the heavy quark momentum diffusion constant 
7 (up to some factors). 

Phenomenological implication of the discretization 
time scale At ~ ^M/CfT is intriguing. Typically, the 
drag force extracted from the experimental data corre¬ 
sponds to Cf 7 ~ [5]. For a charm (bottom) quark in 

the QCP with T ~ 200 MeV, the discretization time scale 
is At ~ 3 fm (5 fm). In the heavy-ion collisions, typical 
lifetime of the QCP is tqgp ~ 10 fm. Therefore the val¬ 
ues of At might indicate that the macroscopic superpo¬ 
sition state remains, at least partially, until the freeze- 
out stage of the heavy-ion collisions. Also the heavy 
quark hadronization time scale 1/Aqcd ~ 1 fm is shorter 
than At and thus in any case the freezeout process may 
be able to resolve the macroscopic superposition state. 
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This would enhance the possibility of recombining heavy 
quark-antiquark pairs into heavy quarkonia and might 
support the statistical hadronization models [ 12 ] when 
heavy quarks are produced in abundance. 

This paper is organized as follows. In Sec. II, I intro¬ 
duce and analyze the basic properties of the heavy quark 
master equation. One can see that Ehrenfest equations 
show kinetic equilibration of heavy quarks. In Sec. Ill, 
I derive the Langevin equation of heavy quarks with 
color. The Langevin equation is not closed in the heavy 
quark phase space: The momentum update depends on 
the heavy quark color state. In Sec. IV, I discuss how 
to interpret the coupling between the dynamics in the 
heavy quark phase space and that in the color space. 
Also the appropriate discretization time scale for classi¬ 
cal Langevin description is discussed. Section V is de¬ 
voted to a summary. Throughout this paper, I adopt the 
natural units, ft. = c = fee = 1, and operators in Hilbert 
and Fock spaces are denoted by bold fonts. 


II. HEAVY QUARKS AS AN OPEN QUANTUM 
SYSTEM 


When one describes a heavy quark in the quark-gluon 
plasma (QGP) as an open quantum system [13], reduced 
density matrix 


Pq( 0 “ T^medPtot (t) (1) 

is a basic dynamical quantity. Here ptot(t) denotes a 
density matrix of a total system of the heavy quark and 
the QGP. The total Hilbert space is composed of a direct 
product of the heavy quark and the QGP Hilbert spaces 
Htot = Hq ® Hmed- By taking the trace over the Hilbert 
space for the QGP (Trmed), one gets the reduced density 
matrix PQ{t) that operates in the heavy quark Hilbert 
space. 

The master equation for the reduced density matrix 
PQ{t) was derived by the influence functional formalism 
[10]. In the leading order of the QCD coupling g and in 
the nonrelativistic limit, ^ it is 

—pQit,x,y) = i ^ pQ{t,x,y) 

+ Fi{x- y)i°-pQ{t, X, y)t - CpFi (6)/5Q(t, 5, y) 

+ h{x-y)-{Vx-Vy)i°‘PQ{t,x,y)i°‘. ( 2 ) 

Here [pQ(t,x,?/)]y = {x,i\pci{t)\y,j) is a reduced density 
matrix in the position and color spaces. The matrix 
is in the fundamental representation of the color SU(A^c) 


^ In this paper, I need not start from the Lindblad form of the 
master equation which preserves the positivity of the reduced 
density matrix. Therefore I neglect all the terms with A{r) in 
Ref. [10]. 


algebra and Cp = (N^ - l)/ 2 Vc. T’i(f) and F 2 (f) are 
given in terms of a function D(f) as 


Fi(r) 

Mr) 



D{f) -H 


v2n(f)\ 
4MT /’ 


-V 


Djr) \ 
AMT) 


( 3 ) 

( 4 ) 


Here D{f) is defined by a thermal two-point function of 
gluons 

^2 ^ oo 


As is clear from the definition, D{f) is an even function 
of f. At ]r| ~ ^IgT, the hard thermal loop resummed 
perturbation theory gives at leading order 


D{f) 



(fk 

(27r)3fc(fc2 + w2)2- 


( 6 ) 


Here wd is the Debye screening mass Wd = {g^T'^l3){Nc + 
Vf/2) for QCD with Nf light flavors. 

From the master equation (2), one can derive time- 
evolution of thermal averaged quantum expectation val¬ 
ues {{0}{t) = Trq {pQ{t)0}, or the Ehrenfest equations, 
such as 



reproducing the consequences of classical Langevin dy¬ 
namics (up to e). The drag force parameter 7 > 0 and 
the correction to equipartition e (jej « 1 ) are given by 

7 = iv"D(r)|,^o, 

^ 4Mr v2D(r)|,^o 

The master equation (2) has several important proper¬ 
ties, which I summarize here. The terms with Fi describe 
thermal fluctuation and that with F 2 describes dissipa¬ 
tion. In the Ehrenfest equations ( 8 ) and (9), the thermal 
fluctuation §^^(1 + e) derives from Fi whereas the 
damping terms —^{{P)) and are from F 2 . 

By taking the trace in the color space, the color averaged 
master equation for pQ{t, x, y) = tr{/5Q(t, x, y)} is derived 
as 

—pQ{t,x,y) = i^^j^pQ{t,x,y) ( 12 ) 

-I- Cp {Fi(x - y) - Fi( 6)) pQ(t,x,y) 

+ CpF 2 {x - y) ■ {VX -Vy)pQ,{t,x,y). 


( 10 ) 

( 11 ) 
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Close to equilibrium, thermal de Broglie wavelength for a 
heavy quark /jb ~ 1 /VMT is much shorter than typical 
correlation length of (electric) gluons Zfluct ~ where 

for |r| > Ifiuct, D{f) ^ 0. Therefore one can take a short 
distance approximation for D{f) ^ Dq + (D 2 / 6 )r^. In 
this approximation, e = 0. By this expansion, one gets 
the Caldeira-Leggett master equation [14]. 

In the following, I derive the Langevin equation it¬ 
self, instead of being satisfied with reproducing the 
consequences of its dynamics. 


III. DERIVATION OF LANGEVIN DYNAMICS 
A. Stochastic master equation 

First, let us investigate more closely how the Fi terms 
in the master equation ( 2 ) are related to the thermal 
fluctuation. Suppose an infinitesimal time step t ^ t + 
dt under a random background ^“(t,x), which rotates a 
heavy quark color state by U^{t,x) = 

The following statistical property of the random field is 
assumed: 

1 (C(t,£)}T = 0. ^ ^ 

Hereafter, I call average over the random field ^ the noise 
average and denote it as (•••)^. Then the density matrix in 
the random field, which I denote as pQ{t,x,y',^), evolves 
as 


PQ{t + dt,x,y-,0 

= U^{t:X)pQ{t,x,f,0 [Udt,y)]’. (14) 

Noting that ^ c< U^{t,x) can be expanded as 

rlf^ 9 

u^(t,x) = i-idtcit,x)i‘^-^[cit,x)i‘^Y+ o(dt^/^) 

yC'F-CiCO). (15) 

In the second line of Eq. (15), I approximate 

dt'^ [^°‘{t,x)F'f ^ = dtCFFi{0) be¬ 

cause the fluctuation around the average is oc dt. Such 
a fluctuation can be neglected, for it does not contribute 
to the master equation after taking the noise average. 
Therefore one can always substitute ^ 

The evolution equation (14) then becomes 

( 16 ) 

= Fi(x - y)i‘^pQ{t,x,y-Oi°' - CFFi(0)pQ(t,x,y-,O 
- i X, y; - pQ{t, x, yiOCit, 2 /)^“] • 

Taking the noise average for Eq. (16) and interpreting 
the noise-averaged density matrix {pQ{t,x,y]Q)^ as the 


original density matrix pQ(t,x,y), one can reproduce 
the thermal fluctuation terms (Fi terms) in the master 
equation (2). Adding the kinetic term and the dissipa¬ 
tion terms {F 2 terms) to Eq. (16), I derive the following 
stochastic master equation: 

Q 

—pQ{t,x,y]0=i^^^j^PQit,x,y-,0 ( 17 ) 

+ Fi{x-y)i°‘pQ{t,x,y;Oi°--CFFi{6)pQ{t,x,y;^) 

+ F2{x-y) ■ (Va: - Vy)t“pQ(t,x,y;Oi“ 

- i[C{t.x)i°‘ Pci{t,x,y-,C) - PQ{t.x,y]OC{t.0''] ■ 

Although I construct the stochastic master equation 
(17) quite intuitively, there also exists a formal deriva¬ 
tion. One can formulate the problem using the influ¬ 
ence functional formalism [15]. Eollowing the notation of 
Ref. [10], the influence functional Sjf can be obtained as 
an expansion Sif = Spot + •S'fluct + Sdiss + •"• The thermal 
fluctuation Sfluct and a part of the dissipation Sdiss can 
be expressed using a stochastic variable equivalent to 
here. See Appendix A for details. 


B. Langevin equation with color 


Using the stochastic master equation (17), one 
can derive stochastic Ehrenfest equations. Since 
Trq {pQ(t;^)(9} depends on the random field it cor¬ 
responds to quantum expectation value in a particular 
realization of which I denote as {0){t;^) or sim¬ 
ply as {0){t) if there is no confusion. The relation 
to the thermal average of quantum expectation value is 
{{0}{t) = ([(( 9 )(t;^)]) 5 . The following stochastic Ehren¬ 
fest equations are derived: 




dt 


{P) 


iP) 

M ’ 


Cf7 

2MT 




(18) 

(19) 


with 


= 7(1 + ^)5{t - ( 20 ) 

In the derivation, I assume that the heavy quark is 
localized as a wave packet. To be explicit, the spatial ex¬ 
tension of the wave packet should be much smaller than 
the correlation length of gluons Zfluct- This is the clas¬ 
sical (point particle) limit for the heavy quark. In such 
a limit, it is natural to assume that the density matrix 
is factorized in the color and the configuration spaces 
Pq{PO ^ Pcolor(t;'C) ® Pconf(t; 0 - ^ In terms of Pcolor, 


^ Although I denote the color space density matrix by an operator 
in the Hilbert space pcoion it is defined only in the Nc x N* 
representation. 
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= Trcoior{Pcoior(i;0*“}- The random force /“ is 
the force in the random potential evaluated at the 
position of wave packet x = {x)(t): 

( 21 ) 

The random force strength (20) follows from the defini¬ 
tions (13) and (21). 

Unlike the conventional Langevin equation, the ones 
derived here in Eqs. (18) and (19) are not in a closed 
form. One needs to know the dynamics of (<“). Using 
the same assumption as made for it is derived 

as 

with 

{Cit)C\t'))T = aS{t-t')6<^K (23) 

Here is the random field evaluated at the position of 
the wave packet C“(l) - = (®)(1))- Note also that 

= 0 which follows from the definitions of 
/“ and C“- The parameter is a = Ui(0) > 0 and is 
structure constant of the color SU(iVc) algebra. Here¬ 
after, 1 write Pcoior(l;C) instead of Pcoior(l;C)- 

The stochastic equation (22) has an interesting prop¬ 
erty: It conserves (t“)(t“). If and only if Pcoior(l;C) is a 
pure state density matrix, (t“)(t“) = | (l - There¬ 
fore Pcoior(l; C) stays a pure state density matrix if it was 
initially. This feature can also be understood concisely 
by reconstructing [/Ocoior(l; C)]ij = (*|Pcoior(l; C)|j) with 

Pcoior(t;C) = ^ + [2(t“)(l)]t“. (24) 

Equation (22) is equivalent to 

^Pcolor(^;C) “ CX^i ^color(^?C)^ “ C^FPcolor(^; C)] 

-^[C“(^)^^/3color(^;C)]• (25) 

This color space master equation is very similar to 
Eq. (16) with substitution x = y = {x){t). Therefore 
Eq. (25) can be derived from random rotation in the 

color space Pcoiorit + dt'X) = 17^(t)/3coior(<; C) [^^cC^)] 
with U({t) = exp[-idtC“(t)t“]. Since the ran¬ 

dom rotation is a unitary evolution, one obtains 
Trcolor ( [Pcolor(^ dtj C)] ) “ Trcolor ( [PcolorC)] ); 

which gives the conservation of (<“)(t“). Note that 
Trcolor 1 [Pcolor(^? C)] | — Trcolor {Pcolor(^; C) } “ ^ SCtS the 
upper limit of Because of the conservation of 

Trcoior I [pcoior(t; C)]^}i T i® natural to demand (only ini¬ 
tially) that Pcoior(i;C) be a pure state density matrix. 

Let us summarize what 1 have obtained so far. Taking 
the classical point particle limit for the heavy quark, 1 


have derived the coupled Langevin equations in the phase 
space and in the color space: 


dr ' M 

A(p) = + r(t)(t“), 

dr ' 2 mt ' ' ^ ' 


(26) 


The statistical properties of the noises are 


'(/r(o/}(0>T = 7(i+e)^(i-o^“%, 

• {C{tX\t'))T = aS{t-t')S<^\ (27) 

,(/r(ocV)>T = o. 


The initial condition for (t“) should satisfy the pure state 
constraint (f“)(t“) = | ^1 - In the above equations, 

the noises /“ and derive from Therefore 1 still use 
the notation of (•••){ for taking the noise average over /“ 
and C“- 

Typical magnitudes of the parameters are 
7=iv^7?(f)|,^o~5"ln(l/ff)T3, 

■ (28) 
AMT V^D{f)\r^o M 

ct = Ei( 0 ) ~ y^T. 

Physically, the drag parameter 7 comes from scatterings 
with both soft and hard momentum exchanges while a 
is dominated by soft scatterings that rotate the heavy 
quark color. 


IV. LANGEVIN DYNAMICS AND 
DECOHERENCE 


A. What is (t“)? 


Having derived the Langevin equation (26), let us in¬ 
vestigate its property in more detail. Hereafter, I regard 
{p) and {x) as classical variables p and x. It is not trivial 
whether to identify (t“) as a classical variable satis¬ 
fying the constraint j^)- Note that is 

now a (IVg - l)-dimensional vector and is different from 
a matrix Actually, this identification contradicts the 
equipartition (9) derived from the master equation be¬ 
cause it results in 


l/iL\ 

dt\2Mj^ Mt(\2m/^ 2^ Cf I' 


(29) 


with rr = |(1“]^)<Cf- Therefore with the identifi¬ 
cation r = (t“) as a classical variable, the heavy quarks 
cannot get thermalized correctly. One easy solution. 
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which turns out to be correct, is to impose a constraint 
= Cy + (<“)(<“), independently of the iden¬ 
tification = (t“). The problem of this quick solution is 
that the connection to Pcoior(i; C) is seemingly lost. I will 
see how to interpret the constraint = Cy- 

To reach a deeper understanding, let us recall the phys¬ 
ical meaning of (t“). By definition, (t“) is a quantum 
expectation value of the heavy quark color charge and 
it couples to the time evolution of classical or macro¬ 
scopic variable p. This indicates that p can be regarded 
as a macroscopic variable in the so-called Schrodinger’s 
cat state. To determine the momentum update, the ac¬ 
tual value for should be the one that is as if 

measured. Otherwise, the noise term only de¬ 

scribes the thermal fluctuation and lacks the quantum 
fluctuation in the color space. 


Wave packets 
in superposition 


Noise field 


A A 

^x 



FIG. 1: Length scales of wave packets in a superposition state 
and the noise field. 


B. “Measurement” of momentum kicks 

In order to take into account the quantum fluctuation, 
I introduce a notation which takes one of the pos¬ 

sible values of the observable O in a measurement with 
respect to a pure state Pcoior(i;C)- Here I examine how 
is measured. First, let us parametrize /“(t) = 
f{t)n°'{t) with normalization n°‘{t)n°‘{t) = 1. To be con¬ 
sistent with the statistical property of /“, the indepen¬ 
dent white noises / and n“ must satisfy {fi{t)fj{t'))T = 
{N^-l)-f{l + e)6{t-t’)Sij and (n“(f)n^(t))T = 

Then r/“(t)f“] takes values as 

L'’ ^ ' Jmeas 

There are Nc - 2 zero modes for in the funda¬ 

mental representation. The probability that a pure state 
Pcoior(i;C) shrinks to each eigenstate depends on the de¬ 
tails of both the pure state and the eigenstate. The quan¬ 
tum fluctuation can also be described using the Wigner 
function. See Appendix B for a brief summary. 

By interpreting the momentum update as 

+ A(0<“]meas: (31) 

one can take into account the quantum fluctuation in the 
color space. In this interpretation, the noise strength due 
to the quantum fluctuation is evaluated as 

/(t)2Treolor{pcolor(l;C) (32) 

Together with the thermal fluctuation, the noise strength 
becomes 

(/(i)'Trcolor {pcolor(i; C) 

= 37(1 + e)d(0)((rt“)) = 3Cf7(1 + e)<5(0), (33) 

and the equipartition of heavy quarks is realized cor¬ 
rectly. In an original derivation, (5(0) = 1/dt with dt being 


a mathematical discretization time scale. The discretiza¬ 
tion time scale dt is shorter than any physical time scales 
of the heavy quark and is thus sent to dt ->• 0. One will 
see, however, that there is an appropriate discretization 
time scale At at which the “measurement” makes sense. 
Therefore, physically, one should take i5(0) = 1/At. 

In the above, the factorization of /“(t) = /(t)n“(t) 
is essential for the simultaneous measurements of 

and In general, /“(t) e -y^^{t,x = {x){t)) 

does not factorize in such a way. However, I need to 
assume the factorization only at the time scale At, dur¬ 
ing which, as one will see shortly, several soft scatterings 
take place. Therefore, the factorization is assumed only 
for the statistical property of the total momentum kick 
accumulated during At. 

C. Decoherence time scale 

Here I show that the quantum fluctuation of /“(t)t“ 
becomes relevant as a consequence of decoherence. For a 
moment, let us get back to the quantum master equation 
for the heavy quark (12). The second line of Eq. (12) de¬ 
scribes decoherence of heavy quark wave function. When 
the heavy quark wave function is localized compared to 
the Debye screening length tfluct ~ l/ffT", Fi{x - y) can 
be Taylor expanded and the decoherence term can be 
approximated as 

^CYj{l + e){x-yfp{t,x,y). (34) 

This term tells us that a wave function of size Ax « 
^fluct loses coherence in a time scale Tdec ~ l/CpdCAx)^. 
If the reduced Planck constant h is recovered, Tdec ~ 

ICyj{P^x)^ indicating the quantum nature of decoher¬ 
ence. 

Let us apply this decoherence time formula to the case 
of the decoherence of the macroscopic superposition in 
the momentum space. The decoherence time gives an es- 
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timate for an appropriate time step At « Tdec in the classi¬ 
cal Langevin dynamics. Figure IV C depicts hierarchy of 
the length scales I consider. When the momentum differ¬ 
ence of the wave packets in the macroscopic superposition 
state is Ap, it propagates to the distance between them 
Ax ~ (AplM)At in time At. The momentum difference 
due to the noise is Ap ~ |/(t)|(2iVc)^^/^At ~ y/UpjAi. 
Therefore I get Ax ~ (At)^/^/M for the dis¬ 

tance between the wave packets. As the time step At 
lengthens, the distance between them Ax increases and 
thus the decoherence time Tdec decreases. The time step 
At of the classical Langevin dynamics should be com¬ 
parable to the decoherence time scale Tdec during which 
the quantum interference between the wave packets is 
destroyed. Equating At « Tdec ~ l/CF 7 (Aa;)^, I obtain 

~~ (ff2iVc)2 ln(l/5)r3 • 


Tcoior ~ 1 /q: ~ which is the interval of soft scat¬ 

terings. Therefore, in a time scale At » Tcoior, PcoiorCO 
is updated, in effect, by a random SU(Ac) matrix. Re¬ 
placing ;Ocoior(i;C) by a random pure state, which subse¬ 
quently shrinks to one of the eigenstates of the 

heavy quark color degrees of freedom are no more dy¬ 
namical degrees of freedom. For a random pure state, 
takes values {l/ 2 , 0 ,-", 0 ,-l/ 2 } with equal 
probabilities. This allows us to treat [»^“(^)<“]i„eas 
a stochastic variable with variance 1/2Ac. In this way, 
the normal Langevin equation without color degrees of 
freedom is derived at time scale of At ~ yjMIC-p^\ 


[ AS = — At, Ap = - pAt -H fj{t)At, 
] M ' ^ 2Mr ^ ’ 

I {Vi{t)Vj{t))T = CftCI + e)6ijlAt. 


(36) 


V. SUMMARY 


With explicit dependence on ft, one can see the quantum 
origin of the time scale At ~ \JhMICpj. In this time 
scale At, the decoherence turns the quantum superpo¬ 
sition in the Schrodinger’s cat state into a probabilistic 
mixture of the classical states. 

There is an alternative view to the time scale At from 
the uncertainty principle [16]. Here let us recover the 
reduced Planck constant ft. Starting from the same 
phase space point, the classical trajectories in superpo¬ 
sition are close to each other by Aa: ~ (Ap/M)At and 
Ap ~ ^CpyAt after time At. Those trajectories get sepa¬ 
rated by AxAp ~ (C'F 7 /M)(At)^ ~ ft at At ~ yfhWJC^ 
and occupy different phase space volume elements. 

The appropriate time step At depends on the mass 
and the momentum diffusion constant of heavy quark. 
^ This is much shorter than the relaxation time of the 
heavy quark momentum At « Treiax ~ MT/Cpy and thus 
is consistent with the above analysis in which the dissi¬ 
pative term with F 2 is ignored. Using this time scale, 
the distance between the heavy quark wave packets in 
superposition is Ax ~ so that Ax « Zfluct 

also holds consistently. Moreover, the wave function is 
much more extended than the thermal de Broglie length 
Aa; » IdB, which is a typical size of a wave packet. 


D. Color space dynamics 

The appropriate discretization time scale for the clas¬ 
sical Langevin dynamics is At. However, the typical 
time scale for the color space dynamics is much shorter 


^ One should also compare At with the hard collision time 
scale Thard ~ l/(tVc - l) 9 '‘ln(l/g)T ~ l/fg^TVcj^T, which pro- 
vides a conventional lower limit of time discretization. Ne¬ 
glecting the logarithms ln(l/p), At > Thard corresponds to 
Yc(M/T)(g2 7V,)2 > 1 , 


In this paper, I studied the Langevin dynamics of 
heavy quarks with internal color degrees of freedom. The 
Langevin dynamics is not obtained in a closed form in the 
phase space but as a coupled dynamics in the phase space 
and in the color space. The coupling causes the macro¬ 
scopic superposition state, or the so-called Schrddinger’s 
cat state, in the updates of heavy quark momentum. The 
classical variable p, which corresponds to the cat, couples 
to the quantum state in the color space. The quantum 
interference in the superposition state is destroyed in a 
time scale At ~ \/ M/Cpy by the decoherence and the 
momentum update should be regarded as one of the pos¬ 
sible outcome by a “measurement”. Note that the macro¬ 
scopic superposition and the decoherence in the Langevin 
dynamics are unique to the non-Abelian interaction. In 
this way, I am naturally led to take At ~ \JMlCp^ for a 
discretization time scale of the Langevin dynamics. It is 
interesting because the physical quantities of the classical 
Langevin equation (M and 7 ) determine the discretiza¬ 
tion time scale At for solving it. At this time scale At, 
the color degrees of freedom are expected to be random¬ 
ized. 

There have been extensive phenomenological studies 
on the heavy quark energy loss and the drag force pa¬ 
rameter ~ has been extracted from the ex¬ 

perimental data [5]. The decoherence of a charm (bot¬ 
tom) quark in the QGP with T ~ 200 MeV proceeds in 
a time scale At ~ 3 fm (5 fm). These values are not 
very small compared to the typical lifetime of the QGP 
in the heavy-ion collisions, tqgp ~ 10 fm. Moreover, it 
is longer than the time scale of hadronization, which I 
roughly estimate to be 1/Aqcd ~ 1 fm. This indicates 
that not only might the classical Langevin equation be 
inappropriate to describe heavy quarks in the fireball, 
but also the freezeout process might be able to resolve 
a heavy quark wave function in the macroscopic super¬ 
position state. Such a wave function is extended about 
Aa: ~ ~ 0.6 fm (0.45 fm) for a charm (bot- 
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tom) quark and may enhance the probability of recom¬ 
bining heavy quarkantiquark pairs into heavy quarkonia. 

To establish a better (semi)classical description for 
heavy quarks in the heavy-ion collisions, it would be de¬ 
sirable to derive a novel kinetic theory for heavy quarks 
in which (i) the time scale below At ~ -y/M/Cpy can 
be resolved and (ii) the superposition of wave packets 
and their decoherence are effectively described. The 
Kadanoff-Baym equation for 4-dimensional Wigner func¬ 
tion can be helpful for this project [17]. 

Finally, let us remark on the similarity between the dis¬ 
cretization time scale At ~ -y/M/Cpy and the formation 
time of induced gluon radiation off an energetic parton 
Tform ~ ~ [18]. Here w and fci are gluon en¬ 

ergy and transverse momentum, and g is a transport coef¬ 
ficient characterizing transverse momentum kicks which 
the energetic parton experiences in the medium. In both 
cases, decoherence determines the time scales. 

For future prospect, it would be an interesting chal¬ 
lenge to simulate the stochastic master equation (17), 
possibly in a simplified version without F 2 terms. With¬ 
out the F 2 terms, the stochastic master equation is equiv¬ 
alent to the stochastic Schrodinger equation [19]. The 
natural time step of the stochastic master equation is 
Tcoior while the classical Langevin dynamics is expected 
to emerge at longer time scale At » Tcoior- Starting from 
a localized wave packet state, one may ask the following 
key questions: 

1. Does the thermal fluctuation with color generate 
a superposition state of wave packets? 

2. How is the interference between the wave packets 
destroyed by the decoherence? 

3. How does pQ{t,x,y;^) acquire the effect of deco¬ 
herence in the corresponding classical phase space 
distribution? 

The answers to these questions will help bridge a gap 
between the quantum and classical descriptions of heavy 
quarks in the high-temperature quark-gluon matter. 
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Appendix A: Thermal fluctuation in the influence 
functional 


In the Markov limit, the influence functional can be ex¬ 
panded as S'lF = fSpot +‘S'fluct +‘S'diss + • • • [10]. The influence 
functional gives a part of path-integral weight for 
the propagation of the reduced density matrix. Explicit 
forms of iSfluct and Sdiss are 






^ ft I 

[ dt f dfxdfyVxD{x - y) 

JtQ J 

■M. «(«)]■/ jM(f] . (A2) 

L J \ / (t jj) 


I 

ir 


Here (p“, j“) is nonrelativistic heavy quark color current 
and labels 1 and 2 denote the Schwinger-Keldysh con¬ 
tour. Integrating by parts, the spatial derivatives in j“, 
'S'diss contains terms of similar structure to 5'fluct- Let us 
call this contribution Then 


"t ^‘^diss 


\ ft y' 

( 


^(Pl) P2)(t,x) 


-1 1 
1 -1 


Pi 

P^2 


(A3) 


Since Ei(f) is positive definite for M » T, the thermal 
fluctuation and a part of the dissipation can 

be expressed using a Gaussian white noise as 


0^5'fluct+^*5^^183 


(A4) 


= exp 




is,y))T = Fi(x-y)S{t-s)S 


ab 


(A5) 


This gives nothing but the stochastic propagation in 
Eq. (14). 


Appendix B: Wigner function and quantum 
fluctuation 

I investigated how to reconcile the kinetic equilibra¬ 
tion of heavy quarks with the Langevin equation (26). I 
found that the quantum fluctuation in the color space is 
essential and proposed an interpretation (31). 

Using the Wigner function, I can also describe the same 
physics and discuss the effect of quantum fluctuation. 
The stochastic master equation (17) with explicit h is ob¬ 
tained by substituting (V^-Vy)/2M ^ ft(V^ - Vy)/2M, 
Fi and ^ The Wigner function is 

defined by 

WQ{t,f,p;0 = J (fse~'’^pQ + • (^ 1 ) 










Using Wq, the color-averaged Wigner function is given 
hy Wq = TrcoiorllUg}. Taking the limit ft ^ 0 in the 
master equation (17), I obtain a Kramers equation for 
Wq: 

-r{t,f) • ^Trcoior{lUQ(t,f,p;Oi“}, (B2) 

with /“(t,r) E -V^“(t,r). The Kramers equation is not 
obtained as a closed form of IUq but rather depends on 
the color state of IUq. Let us make the equivalent as¬ 
sumptions which have been made previously: 

1. IUq is localized in the phase space. This allows us 


to replace /“(ftr) by /“(t). 

2. lUq can be factorized into the Wigner function 
in the phase space and that in the color space: 
Wait,r,p-0 = Wait, (t ; 0. 

3. The momentum kick at the position of the heavy 
quark can also be factorized: /“(t) = f{t)n°‘{t). 

In order to take the quantum fluctuation into account, 
the last term of the Kramers equation can be understood 
as 

(B3) 

where is evaluated by the color state of 

The dynamics of lUQ°'°‘'(t; is the same 
with that of /3coior(i;C) in Eq. (25). 
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